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Abstract. 

We theoretically investigate cross-talk in hyperfine gate control of donor-qubit 
quantum computer architectures, in particular the Kane proposal. By numerically 
solving the Poisson and Schrodinger equations for the gated donor system, we calculate 
the change in hyperfine coupling and thus the error in spin-rotation for the donor 
nuclear-electron spin system, as the gate-donor distance is varied. We thus determine 
the effect of cross-talk - the inadvertent effect on non-target neighbouring qubits - 
which occurs due to closeness of the control gates (20-30nm). The use of compensation 
protocols is investigated, whereby the extent of crosstalk is limited by the application 
of compensation bias to a series of gates. In light of these factors the architectural 
implications are then considered. 
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1. Introduction 

There is a growing interest in solid state proposals for quantum computation, owing 
to their promise of a scalable technology given the connection with existing nano- 
fabrication techniques. An important scheme in silicon was proposed by Kane pQ, which 
although difficult to realise, has a number of advantages, including long decoherence 
times of both 31 P nuclear and electron spins [2J. Dopant based proposals, in general 
propose qubits to be encoded via the spin state of the donor-nucleus or electron 
|3j or donor-charge state [I] of a phosphorus atom in a silicon substrate. Single qubit 
control in the Kane case is achieved by the application of gate potentials (derived from 
surface electrodes), which shift the hyperfine interaction of a given target qubit in or 
out of resonance with a global transverse AC magnetic field. Two-qubit operations are 
performed by drawing together the two donor-electron wave functions, such that the 
exchange interaction is modified 0EIIII- This requires the donors to be separated at 
distances of 20-30nm. 
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The standard paradigm for quantum computing requires the ability to individually 
control single qubits in a register of many qubits, as well as to implement controlled inter- 
qubit interactions. The proximity of neighbouring qubits results in poor localisation of 
gate potential, giving rise to a cross-talk effect on non-target qubits. This cross-talk in 
quantum control renders selective addressability for single qubit control problematic. A 
similar problem arises in other dopant based proposals [13 El E] , superconducting qubits 
and in quantum dots |10t lllj. Analogous to this control problem, measurement cross- 
talk in such architectures may also be a problem (similar to the SQUID equivalent j!2j). 
Although both cross-talk and decoherence processes cause errors in qubit operations, we 
wish to distinguish the two. Decoherence is a fundamentally non-unitary process caused 
by interaction with the environment, whereas cross-talk leads to unwanted unitary 
evolution of non-target qubits during an operation. Given the stringent requirements 
of control precision that fault-tolerant quantum computation imposes at the 10 -4 
level [13J, determination and elimination of cross-talk is a key issue for scalable quantum 
computation. 

This paper is organised as follows: In Sec. 2 the physical architecture and 
considerations for modelling is described. In Sec. 3 the calculation of donor-electron wave 
functions and hyperfine coupling is presented with an emphasis to illustrate the cross- 
talk problem and the fidelity of single qubit addressing. In Sec. 4 the gate compensation 
to improve selectivity is described and considers alternative configuration of. 

2. Modelling the Kane architecture 

A schematic of the system we consider is shown in Fig. ^ This architecture consists 
of 31 P (spin I = \) atoms doped in purified 28 Si (1 = 0). The dopants are separated 
by 30nm and embedded at a depth of 20nm below the 5nm oxide layer. A series of 
surface gates referred to as A (aligned above the donor) and J (in between the donors) 
are placed on a 5nm thick oxide layer. The gates are lOnm wide and separated by 5nm. 
A backgate is connected to the substrate by a graded n-doped layer such that it forms 
an ohmic contact withO the substrate. 

Each 31 P atom forms four covalent bonds with neighbouring 28 Si atoms, with the 
fifth valence electron loosely bound to the P + ion, with a Bohr radius of a ~ 2nm. 
Operations to be performed on the system are governed by the effective spin Hamiltonian 
for the 2-donor system, given by: 

H*t** = H z + H A + H J + H aB , (1) 

where Hz is the Zeeman term due to a constant magnetic field applied to the entire 
system 

H z = ^ B B{a\ e + a 2e ) - g nf i n B(al n + of), (2) 

Ha is the contribution of the hyperfine interaction that couples nuclear and electron 
spins 

H A = A x o Xe ■ a ln + A 2 a 2e ■ a 2n , (3) 
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Figure 1. a) A Kane architecture based model, b) schematic configuration of gates 
and donors. 



Hj is the contribution of the exchange interaction coupling neighbouring electrons 

Hj = J 12 a le ■ cr 2e , (4) 

and if ac is the contribution of the transverse rotating magnetic field to the Hamiltonian 

#ac = -9niinBJia^ + al n ) cos(wt) + (<rj n + of ) sm(wt)\ (5) 
+ n B B &c [(al e + a 2 x e ) cos(wt) + (a l y e + af) sin(wt)]. 

In the above expressions fiB is the Bohr magneton, /i n is the nuclear magneton, g n is 
the Lande factor for 3 1P, B is the magnitude of the constant magnetic field, B ac is the 
magnitude of the global transverse AC magnetic field, w the frequency of the transverse 
field and a x ,a y ,a z the Pauli matrices operating on electron and nuclear spins. A and J 
stand for hyperfine and exchange couplings respectively. For a detailed description of 
the Kane architecture refer to [Tl^ I15j. 

The coupling parameters A and J are controlled via A and J gates, in order to 
perform qubit operations. The problem of cross talk can be summarised as follows; 
rather than being dependent on a single control gate, the Hamiltonian parameters A, J 
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are functions of biases applied to all the control gates in the vicinity: 

A = A (V bias ). (6) 

where Vbi as = {Va , Vj h Vai, •••} is the set of voltage biases applied on the gates (Fig.Q 
b)). This dependence and the inherent crosstalk has not received investigative attention 
so far. 

3. Gate control of donor-electron wave functions 

The couplings A and J depend on the donor electron wave function, which is 
determined by the solution of the electronic wave function Hamiltonian. There have 
been a number of approaches to the calculation of the gate induced shift in the 
hyperfme coupling strength at various levels of sophistication of both electron wave 
functions and description of control field, including: perturbation theory with envelope 
hydrogenic wave function for circular and strip electrodes |lfij . envelope hydrogenic basis 
diagonalisation using numerically determined nano-gate potentials |17| ITS], variational 
effective mass theory in uniform fields [TH]. tight-binding calculations with uniform field 
PU] . self-consistent calculation of 3D Poisson and DFT taking into account the six 
valleys of Si band structure [21] and numerical expansion in a basis of silicon Bloch 
functions in a uniform field f2ZJ. Despite differences in detail, there is good agreement 
between these various methods on the hyperfme response to external electric fields. A 
hydrogenic approach has proven to be sufficient for the calculation of the hyperfme shift, 
whereas the details of the potential field generated by the control gate are vital to the 
description of the cross-talk problem. Therefore, in this first determination of the extent 
of, and possible correction procedures for cross-talk in the quantum control problem, we 
adopt a hydrogenic envelope approximation for the donor electron wave function, but 
determine the electrostatic fields numerically to get a reliable description. 

In this analysis we ignore the effects of the image charges of donor nuclei and 
electrons in the silicon-oxide barrier, as well as the effect of the magnetic fields on the 
electron wave functions. Both of these issues are discussed in |2S], and will alter the 
energies of the bound-states of the inerface well that is formed beneath the control gate. 
This will change the voltage at which it becomes energetically favourable to occupy 
these surface states for all donors, and will not produce a relative energy shift between 
the donors. We also ignore any inter-donor coupling, to include this would complicate 
the calculation considerably as it would be necessary to calculate the ground-state of a 
many electron system. We expect that inter-donor coupling will only play a role when 
the electron wavefunctions have been significantly deformed, and we are not interested 
in any regime in wich more than one donor electron wavefunction is distorted. 

The ground-state single-donor wave functions are calculated by expanding in the 
basis of scaled hydrogen-like orbitals. The single-donor electronic Hamiltonian in the 
presence of an external potential generated by a bias profile Vbi as applied to a surface 
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gates, is given by 



H „ Q e 2 



^eiect = _ v 2 - — + U(r, V bias ) + C/ bairier (r). (7) 

2m 47reoe S jr 

In the basis of scaled Hydrogen-like orbitals we compute the matrix elements 

(n,l,m\H elect \n',l',m') = J K,i,Jf)[W> V bias) + U hailie r(?)\*n>,i>, m >(f f )dr + 

+ S n nf5ni8 mm /—^- , (8) 

where n, I, m are the electronic, orbital and magnetic quantum numbers and E n = Ei/n 2 
is the energy of the orbital. In the scaled donor units, as = 2nm and E\ = — 45meV. 
The potential U(r, Vbi as ) is composed of the externally applied electrostatic potential, 
and t/barrier(r) * s the potential barrier (3.25eV) across the semiconductor-oxide interface, 
which arises due to the relative energies of the conduction bands of the two materials. 
The potentials in the device are calculated by solving the Poisson equation using ISE- 
TCAD [23J. The calculations are performed for a device temperature of IK, where 
convergence is obtained by bootstrapping room temperature data carefully down to 
required temperature. The matrix elements are calculated numerically using the Monte 
Carlo integration, and the Hamiltonian is diagonalised to find the donor-electron ground- 
state as a function of the applied gate biases. 

3.1. Fidelity of single- qubit addressing 

Single qubit operations are performed by controlling the donor nuclear-electron hyperfine 
coupling A, via gate potentials. Both A, and the exchange coupling J, depend on 
the wave functions of the donor-electrons. Application of a surface gate bias deforms 
the donor electron wave functions, thus altering the exchange and hyperfine coupling 
parameters. The hyperfine coupling between the nucleus and the electron depends on 
the overlap of their respective wave functions. The probability density of the nucleus is 
treated as a delta function, so the coupling strength A is given by 

^4(V bias ) = y/i B £n/wU^(0,V bias )| 2 , (9) 

where /i s ; is the permeability of silicon and \l/(0, V b i as ) is the donor electron wave function 
evaluated at the nucleus [TS]. Due to complications arising from the core electrons of 
the donor, it is difficult to calculate A directly, rather we calculate the relative change 
in coupling, as a function of applied field [22] : 

*(V«-) = l ^' { o,of A{0) ' (10) 



where A(0) is determined from experimental data 

Single qubit operations are implemented by bringing the target qubit into resonance 
with an oscillating transverse magnetic field. Changing A changes the resonance 
frequency of the donor nucleus. To do this selectively requires the ability to sufficiently 
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alter the hyperfine coupling strength of an individual donor, while keeping the other 
donors sufficiently far off resonance such that their spin state remains unaltered. 
As discussed, the close spacing of the qubits makes unique addressing problematic; 
individual hyperfine control is hindered by the gate potential cross-talk as shown in 
Fig. El a), where we have plotted the hyperfine coupling coefficient as a function of gate 
bias for several bits. 



a) b) 
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Figure 2. a)Response of the hyperfine coupling of an array of donors (as in Fig.H]b)), 
b) Comparison of the Hyperfine coupling constant 'A' of target donor Qq at various 
depths. 

The results in figure Efb) are in good agreement with previous results [El HU GDI EH] 
and show the donor-electron response when placed extremely close to the oxide layer. 
The strength of the hyperfine coupling is also affected by donor depth, due to effects of 
both gate potential and oxide barrier. While the hyperfine coupling at donor depths > 
12nm from the oxide follow a regular pattern influenced by gate potential, very shallow 
depths show pronounced effects due to the oxide layer. At very shallow depths, we 
find that the hyperfine coupling initially increases with the applied positive bias before 
decreasing. This is due to the fact that when the donor is close to the oxide, the donor- 
electron at zero bias is pushed away from oxide, a positive bias restores A to the normal 
value and then A decreases as the wave function is pulled further towards the oxide. 

As the gate voltage is increased, we find that the electrons from all donors are 
perturbed, as shown in Fig. El It is found that for positive applied biases the electron 
wave functions of neighbouring donors are perturbed at lower voltages than the target 
qubit, while for negative biases the behaviour is more intuitive, with the target electron 
perturbed at a lower bias. 

This can be explained as follows; the positive gate potential pulls the electron wave 
functions towards the oxide, while the oxide restricts the wave functions from entering 
it, forming a potential well at the silicon-oxide interface, underneath the gate. The 
depth of this potential well is the same for all donors, however donors nearer to the gate 
have the potential well of their ionic core lowered by the application of the field, more 
than for donors further from the gate. Therefore, it becomes energetically favourable 
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for the electrons of these distant donors donors to occupy this interface well at a lower 
potential than is the case for donors closer to the gate (Fig. 0J). It should be noted 
however, that the time scales on which these electrons tunnel into this interface well 
increases as the donors become further from the gate, as discussed in ref |25J. Therefore, 
it may be possible to use a kind of temporal selection whereby the gate pulses are timed 
such that they allow adiabatic evolution of the target donor electron into the interface 
well, but are too fast for the neighbouring donor electrons to tunnel. The design of such 
pulses with the requisite fidelity for quantum computation would be a formidable task, 
involving solution of he time-dependent Schrodinger equation in a very large basis of 
states. Here we restrict our attention to designing "ground-state" protocols, in which we 
assume that the pulses are applied in such a fashion that each donor electron is allowed 
to evolve adiabatically in its ground electronic state. 

In the case of a negative applied potential, the potential wells of the donor cores are 
raised in energy by the gate potential, the donors close to the gate are affected more by 
this than those further away, and so are perturbed at a lower gate bias. However, this 
results in a very weakly bound donor electron state, which may be subject to significantly 
increased spin-orbit coupling, leading to increased dephasing rates. Additionally, such 
weakly bound states are more susceptible to non-reversible ionisation. We are, therefore, 
hesitant to propose this as a practical method for implementing gate operations. 

3.2. Effective gate errors 

Single-qubit operations are implemented by bringing the energy difference between spin 
states of the nucleus into resonance with a globally applied, oscillating transverse field, 
for a time t, which determines the angle of the rotation. The resonance frequency u res , 
of the nucleus is given by[T]: 



At time t, the probability of measuring the i nucleus in its spin flipped state is 
given by Rabi's formula [2*7j : 



where v is the frequency of the field, z/* es is the resonant frequency of the nucleus and 
7 = 2g n fj, n B ac /h. For the target qubit Q , we define the error as e = 1 — max[P (t)], 
which gives the error due to under- rotation of the qubit. In the case of the non-target 
qubit (off-resonant with the RF field in general), we define the error as e^o = m ax[Pj(t)], 
which quantifies any unwanted rotation of the qubit. In what follows we consider the 
error in implementing spin flip on a particular target qubit, we assume that there is no 
error associated with the timing of the pulse and so set the sine function in the above 
expression to unity. 



hu ies = g Q /j n B - fx B B + yj AA 2 - (g n fi Q B + n B B) 2 + 2 A. 



(11) 




(12) 
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Figure 3. Donor-electron probability density under perturbation when +1.0V is 
applied on A gate (all other gates are grounded) for a) target donor Q and c) 
neighbouring donor Q±. Negative gate voltages produces different behaviour as shown 
when A gate is biased at -0.8V b) for Q and d) for Q%. 



Studies in fault tolerant quantum computation provide a rough estimate of the 
allowable error in any gate operation, which is no more than one part in 10 4 In 
the case of single qubit rotations, the error is composed both of contributions associated 
with the target qubit being out of resonance, as well as neighbouring qubits being 
insufficiently off resonance so that they undergo a partial rotation. 

The hyperfme coupling is least sensitive to gate voltage fluctuations when 
dA/dV=0, which is the case when either A = A(0), or, with the exception of very 
shallow depths, for A = [28J, as seen in (Fig. |2]b)). This suggests that the frequency 
of the global field should be chosen to be the nuclear resonance frequency of the target 
qubit in one of those cases. We set the AC field to be resonant with the A = case, 
in the first instance. Plotted in Fig. a) is the gate error for each qubit as we change 
the target gate bias, at this frequency. We find that the required error rates are not 
achieved for positive gate biases. 

To determine optimum frequency for the global field, for the implementation of 
single qubit operations, we allow the global field to be the nuclear frequency of the 
target qubit [y = z/ res ) and sweep the gate bias for the relevant range. These results 
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a) b) 




Figure 4. a) The effect of gate potential increases with gate-donor distances, b) 
Schematic (not to scale) illustration of gate potential wells: i) potential well created 
by donor-electron in the absence of gate potential, ii) lowering depth of donor- electron 
potential well due to gate voltage, iii) lowering of depth is less at increased distance 
from gate, hence S\ > Sq. 



a) b) 
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Figure 5. Evaluating performance of single qubit operation on a qubit array : a) 
when resonance frequency is chosen to be nuclear frequency of target when A$=§, b) 
when resonance frequency is chosen to be the nuclear frequency of the target donor 
(Qo does not have any error and hence not plotted). 

are plotted in (Fig. 0b)), and show that positive biases do not give a practical window 
of operation. The spike at around IV is due to the crossover of the hyperfine coupling 
between target and neighbouring donor-electrons (Fig. |2]a)). It is clear that some kind 
of compensation is required to reduce the cross-talk. 
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4. Compensation for cross-talk induced gate errors 



The diffuse nature of the electric fields produced inside the device by the application 
of a bias voltage to any particular control gate makes selective addressability of qubits 
problematic. In this section we investigate a simple linear compensation scheme to 
address this issue. Through the application of compensation biases to neighbouring 
gates we minimise the electric potential in the vicinity of non-target donors whilst 
maintaining as much as possible in the vicinity of the target qubit. Therefore we design 
a bias profile such that both the potential, and x-component of the electric field, at the 
positions of the non-target qubits be zero. The required compensation gate profile is 
achieved by the following procedure. Assuming linearity of potentials in the device, the 
potential and it's x-gradient, at the donor sites can be described as: 

n = MV bias (13) 

where V bias refer to the gate biases, and II refer to potentials and their gradients at the 
respective qubits: 

U=(U Qo ,U Qo ,U Ql ,U Qi ,...) T . (14) 

Here 

Uq = U(r Qo ;V hias ) = u Ao (f Qo )V Ao + uj 01 (f Qo )Vj 01 + ... , (15) 

and 



UA (r Qo ) = U(r;{V Ao = 1,0,0,...}) 



(16) 



with u Ao (rQ ), the potential at qubit Q due a bias V Ao = IV. The matrix M is 
constructed by solution of the Poisson equation using a commercial software package 
TCAD, such that Eqn. ED can be written as follows 

( u A (r Qo ) u Joi (r Qo ) u Al (r Qo ) ■■■ u An (r Qo ) \ ( 
UA (r Ql ) 

u' Ao {r Ql ) 



( U Q 

Uq 

U'n 



\ 



u 
w 



Uj 01 (r Ql ) 



\ 



u' A 



V Ao 
Vj 01 

v Al 



V 



(n — l)n 

V An 



UA {r Qn ) 
u' Ao (r Qn ) 

where U' is the derivative of potential with respect to x-direction. We require a value 
of Vbias, such that gives a II with Uq = 1, and all other elements zero. The gradient 
of the potential at the target donor does not need compensation due to symmetry, and 
hence is ignored. This is obtained by matrix inversion; 



V 



bias 



M~ l U. 



(18) 



Applying the above compensated bias profile produces a potential that is more 
localised around the target donor (Fig. EDb)). The compensation assumed linearity 
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Figure 6. Potential profiles obtained from Poisson solver: a) one gate is biased while 
the rest are grounded b) after compensation, the potentials directed at the target is 
more localised. 



of the medium, which would be the case in the absence of charge carriers. However 
non-linearities in potentials introduced by the presence of P-dopants (10 11 /cm 3 ) in the 
semiconductor, result in slight deviation from ideal compensation at all donor-sites, as 
can be seen in (Fig. 0a)). 



a) b) 
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Figure 7. a) Potential profile before and after application of compensation biases to 
other gates, b)Hyperfine coupling of donor array after compensation with respect to 
V Ao - 

The response of electron-nuclear hyperfine coupling to compensation biases is 
plotted in Fig. 0b). One of the consequences of compensation is that larger voltages are 
required to affect the donors, with a 2V bias required to produce the same effect as IV 
in the uncompensated case. This is due to the application of counter-acting voltages on 
the neigbouring gates, which have the effect of reducing the overall potential within the 
device. In (Fig. |BJ) we show the probability density for target and neighbouring donor 
electrons with a compensated bias of 2.6V applied to the target gate Va . This shows 
that, in contrast to the uncompensated case, the target qubit feels a greater effect than 
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Figure 8. Wave function response of a) Qo, b) Q\ after compensation. 
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Figure 9. Error probability : When resonance frequency is chosen to be nuclear 
frequency of target when ^4o=0 , b)When resonance frequency is chosen to be the 
nuclear frequency of the target donor (Qo does not have any error by definition in this 
case and hence not plotted). 



Evaluating the fidelity of a single qubit rotation on a compensated qubit array, 
shows a marked difference to the uncompensated case (Fig. EJ). In the positive voltage 
regime, there exists a small window, at about 2V - 2.5V, in which the required fidelity 
is achieved (Fig. El b)). However, this may still be too constrained. 

In light of this we consider an architecture which consists of a series of well 
separated two-qubit cells, where quantum information can be transported between cells 
via coherent electron transport [HI EI], as sketched in (Fig. [TUla)). We have calculated 
the fidelity of single qubit operations for both the uncompensated, and a compensated 
gate profile, which we show in (Fig. Elb)). This gives a more robust compensated 
voltage regime in which to perform single qubit operations. 
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Figure 10. a) Donors organised as two-qubit cells interacting via quantum transport 
mechanisms, b) Error probabilities of rotating non-target qubit when the transverse 
RF is fixed at the resonant frequency of target nucleus. 



5. Conclusion 

Practical quantum devices require scalable architectures where crosstalk would have to 
be minimised. The feasibility of the conventional Kane architecture has been analysed in 
light of such problems that arise due to close proximity of qubits. The concept of cross- 
talk in dopant-qubit architectures has been introduced and the physics behind this effect 
is explained. We have shown that cross-talk can be reduced using gate voltages based 
on linear compensation techniques, however, the use of more sophisticated quantum 
control techniques may be also be beneficial. 

In response to the various issues for scale-up of the original Kane proposal, which 
includes the problem of cross-talk analysed here, an alternative quasi 2D donor based 
architecture has been proposed jHOj- The initial analysis of cross-talk for the two donor 
interaction cells shows that the problem can be more effectively compensated in such 
an architecture. 
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